Impurity correlations in dilute Kondo alloys 
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The single impurity Kondo model is often used to describe metals with dilute concentrations 
(m) of magnetic impurities. Here we examine how dilute the impurities must be for this to be valid 
by developing a virial expansion in impurity density. The 0(nf ) term is determined from results on 
the 2-impurity Kondo problem by averaging over the RKKY coupling. The non-trivial fixed point 
of the 2-impurity problem could produce novel singularities in the heat capacity of dilute alloys at 
0{n 2 t ). 
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I. INTRODUCTION 



The physics of the single-impurity Kondo model, the 
problem describing a magnetic moment interacting with 
a sea of conduction electrons, is well understood 13]. At 
high temperatures the moment is free. The spin sus- 
ceptibility of the impurity obeys a simple Curie- Weiss 
law, and there is a large additional term in the entropy, 
S = fcsln2. As the temperature decreases, the impu- 
rity spin disappears as it forms a singlet with conduction 
electrons. The characteristic energy scale at which this 
happens is the Kondo temperature, which is exponen- 
tially small at small couplings, Tr ~ EFexp(—l/(Jp)), 
where J is the s — d exchange constant, p is the den- 
sity of states. The 1/T divergence of the spin suscep- 
tibility is then effectively cut off by this scale at low 
temperatures. The singlet state formed by conduction 
electrons and impurity can be described at low temper- 
atures T <C Tk as a Fermi liquid with enhanced density 
of states, bv ~ tii/Tr, where rii is the concentration of 
impurities. 

The single-impurity Kondo model is often used as a 
way to obtain insights into the nature of the ground state 
for the Kondo lattice problem, a toy model for a num- 
ber of rare earth compounds, such as heavy fermions. 
Enhancement of the density of states due to the Kondo 
physics explains qualitatively the heavy masses in the 
heavy fermion materials. The local moments in a Kondo 
lattice interact through the electron-mediated RKKY in- 
teraction, which competes with the Kondo physics. The 
low temperature behavior of the Kondo lattice is often 
thought to be determined by this competition. Conven- 
tional mean field models of the heavy fermions , how- 
ever, often neglect the RKKY interaction, since it only 



appears as a correction to the mean field behavior, and 
therefore is difficult to calculate. 

The simplest model that captures the physics of this 
competition is the two-impurity Kondo model UK . De- 
pending on the relation between the RKKY interaction 
and the Kondo temperature and the presence of particle- 
hole symmetry breaking, the impurity spins can each get 
get compensated by conduction electrons, form a singlet 
with each other, or a combination of both. A surface of 
Fermi liquid fixed points occurs in the problem as well 
as an unstable non-trivial critical point associated with 
critical behavior of the thermodynamic functions [p|dTl]| . 

Despite extensive theoretical work on the 2-impurity 
problem, we are not aware of any effort to apply the 
results in a systematic way to experiments. The appro- 
priate way to do this, for dilute alloys, was pointed out 
in 1970 by Larkin and collaborators |12|. One should 
develop an expansion of thermodynamic quantities in 
powers of the impurity concentration, rii, a virial ex- 
pansion. The pioneering work on this subject Jl2|| only 
considered the high-T limit where the Kondo interaction 
could be ignored and only the RKKY interaction consid- 
ered. At T of 0(Tk) or lower the Kondo interaction must 
be included. The 0(nf) term in the virial expansion is 
completely determined by properties of the 2-impurity 
problem appropriately averaged over the strength of the 
RKKY interaction. This calculation gives an estimate 
of the characteristic impurity concentration at which the 
single impurity Kondo model breaks down. Analysis of 
the A-impurity models with N > 2 could provide a fur- 
ther insight into the structure of this expansion. 

We note that deviations from linear dependence of 
thermodynamic quantities on impurity concentration do 
not arise solely from the mechanism considered here. The 
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Kondo temperature depends exponentially on the ex- 
change integral, and changes in the composition of an 
alloy can alter it significantly. Such corrections, how- 
ever, can be incorporated in the single-impurity picture, 
even though they alter the linear dependence of the ther- 
modynamic and magnetic properties of a dilute alloy on 
concentration. 

In the next section we discuss the subtle issue of length 
scales in the Kondo problem. Section III sets up the 
virial expansion, following £[§]. Section IV reviews some 
relevant results on the 2-impurity Kondo problem. The 
0(nf) term is then considered for T ^> Tk in Section V 
and at T < Tk in Section VI. Quantitative results would 
certainly require accurate numerical solution of the 2- 
impurity problem for all values or the RKKY coupling. 
Here we just apply the large-N approximation and make 
some qualitative remarks on the effect of the non-trivial 
critical point. Section VII contains conclusions. 

II. LENGTH SCALES IN THE KONDO 
PROBLEM 

There are (at least) 3 different ways of estimating a 
lower bound on the average separation of impurities nec- 
cessary for the single impurity model to be valid. One 
could require that each impurity be far from the others 
compared to the size £k of its screening cloud. Alterna- 
tivelyone could require that the density of electron states 
within an energy Tk of the Fermi surface be at least as 
large as ni (the Nozieres exhaustion principle). This ap- 
pears to be a neccessary (although not obviously suffi- 
cient) condition for the screening cloud wave-functions 
from each impurity to form an approximately orthogo- 
nal set. Or, finally, one could require that the average 
RKKY interaction be small compared to Tk- 

As a mechanism of Kondo screening, the picture of spin 
exhaustion cloud of exponentially large size £k ~ vf/Tk 
(vf is the Fermi velocity) has often been adapted 0. 
From the RG point of view, clearly, this length scale is 
present in the problem, since the low-energy Fermi liq- 
uid theory is only valid for k — kp <C Tk/vf, or at 



So far we have only used the spherical symmetry of 
the screening cloud wave-function and hence its Fourier 
transform. Now we assume: 

\^(k)\ 2 ^(U/k 2 F)f[(k-k F )U}, (2) 

where the scaling function f(y) obeys the normalization 
condition: 

J dyf(y) = 2tt 2 , (3) 



distances r 3> £,k away from the impurity. Direct calcu- 
lations of the screening cloud profile |H for the spin-spin 
correlators of the single-impurity Kondo model demon- 
strate that the screening cloud of conduction electrons 
indeed forms at the distance scale ~ £k. Note that this 
is a more dynamical type of a screening than that which 
occurs for charge impurities in a Fermi liquid since it 
involves a linear combination of states where the impu- 
rity spin and the screening electron spin are in cither 
an up-down or down-up configuration. In particular, the 
finiteness of the susceptibility at T — > should not be at- 
tributed to a static conduction electron polarization can- 
celling the impurity spin polarization. Rather it results 
from the tendency of the impurity to form a singlet with 
the screening electron. This tendency is very well illus- 
trated by a calculation of an equal-time correlator, which 
provides an instant snapshot of the system, and the zero- 
frequency spin correlator, or the spin susceptibility [fl4|| . 
Indeed, it has been shown that the magnetic moment of 
conduction electrons exactly compensates that of the im- 
purity at the same moment in time. Yet, there is no net 
polarization of the conduction electrons (as illustrated by 
a calculation of zero- frequency spin susceptibility), the 
fact also known as the Anderson-Clogston theorem. 

Existence of an exponentially large length scale could 
have potentially strong effect on the theory of alloys with 
magnetic impurities. Indeed, typical Tk ~ 10 K and 
Ep ~ lOeV makes £k ~ 10000a, where a is the lattice 
spacing, much larger than typical distance between two 
impurities. This issue was addressed in one dimension for 
Luttinger liquids with magnetic impurities p5[ , where it 
was found that the crossover happens at 7ij ~ 1/Ck- At 
higher dimension, however, no such crossover was found 
experimentally. The physical reason for this is that the 
screening cloud wave function has an oscillatory char- 
acter. Let us make a simple estimate of the overlap of 
the screening clouds from two separate impurities, with 
one impurity located at a distance R apart from the 
other. We shall make no particular assumptions about 
the screening cloud wave function except that it is made 
up of Fourier modes with \k — kp\ less than or of order 
of 1/^k- Thus the overlap is: 

(1) 



in order that O(0) = 1. Thus we get: 

O(R) = (l/2ir 2 ) J dyf(y) sm[k F R + {R/£ K )y]/k F R. 

(4) 

For R<^ £k this reduces to: 



0(R) = I d 6 rr {\r + R/2\)i,{r-R/2\)= I -^ e *^(fc)| 2 = [°° dkk 2 \^(k)f MkR) 



(2tt) 3 Ir ' n J Irv n 2w 2 kR 
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O(R) = sm(k F R)/k F R. 



(5) 



fi = Fi 



Thus the overlap is small for k F R^> 1. This calculation 
can be easily extended to the 1 or 2 D case. For 1 D the 
overlap is essentially cos(2k F R) for R -C £,k- It doesn't 
get small until R £k- In 2 D it is a Bessel function, 
J (k F R) « ^2/irk F Rcos(k F R - tt/4) (for fc F i? > 1), 
again becoming small for k F R 1. We see that in di- 
mensions higher than one the overlap of the two screening 
cloud wave functions is suppressed by l/y/k F R(2D) or 
l/k F R(3D), small factors at large enough inter-impurity 
separation. In ID this overlap is 0(1). Thus the dimen- 
sionality of the problem could be the main reason why 
the large Kondo scale is absent in alloys. 

The well-known Nozieres' exhaustion principle states 
that for each impurity spin there should be a conduc- 
tion electron in the vicinity ~ Tk of the Fermi energy 
that screens it. This produces the following estimate for 
the inter-impurity separation: R ~ (Crt&f ) 1/ ' 3 in 3D, 

R ~ \JiKkp 1 in 2D, and R ~ £ K in ID. 

On the other hand, one can argue that the single- 
impurity model stops working when the RKKY interac- 
tion at the average inter-impurity distance becomes com- 
parable to the Kondo temperature, Jrkky ~ Tr. This 
gives a somewhat higher estimate for the concentration 
at which one needs to account multi-impurity effects (or 
lower estimate for the average inter- impurity distance), 

R ~ (\ 2 Z K kp 2 ) 1/3 in 3D, R ~ \j \ 2 Uk F 1 in 2D, and 

R ~ \ 2 £,k in ID, because of its dependence on the small 
coupling constant, A = (Jp). 



III. THE VIRIAL EXPANSION. 

We start with the ordinary s-d Hamiltonian for a col- 
lection of magnetic impurities in an electron gas located 
at arbitrary positions in space: 

H °<* = ~Wn E ' S LXe l(k - k '^ (6) 



2N 

The free energy for N impurities is determined by the 
formula: 

F 1)2j ... )W = -Tln(exp{- / ' H sd (r)dT}) (7) 

Jo 

To obtain a virial expansion, one can introduce the quan- 
tities / defined by recurrence relation: 



r---m "i" fij"-ni (p) 



where the summation is carried out over different sets of 
indices ij ■ ■ ■ n. The function / vanishes if the distance 
between any two impurities tends to infinity. For exam- 
ple, 



fij — Fij (Fi + Fj). 



(9) 



Averaging Eq.(g) over the distribution of the impurities 
and going to the thermodynamic limit N — > oo, we obtain 
an expansion of the free energy in powers of the density. 
For magnetic impurities interacting with an RKKY in- 
teraction such expansion was carried out in Ref Q . The 
nf term in the Kondo regime without the RKKY term 
was calculated in Ref. , although the influence of the 
RKKY interaction on this result was discussed in a later 
paper. The ordinary Kondo term is proportional to the 
density of impurities, n,; = Ni/V. The n? contribution 
takes the following form: 



F {2) = N, 



dK(F 2i (R) - 2F, 



Kondo ) 



(10) 



Here F 2 i(R) is the free energy of the two impurities sepa- 
rated by a distance R, FKondo is the usual single- impurity 
Kondo term. Thus the n 2 correction is competely deter- 
mined by the two-impurity physics. Eq. (|To|) and similar 
expressions for thermodynamic functions (derivatives of 
F) could be useful for numerical determination of the 
temperature dependence of the nf term in thermody- 
namic functions from the results on the two-impurity 
model. 



IV. THE TWO-IMPURITY KONDO MODEL. 



Before calculating the nf term let us discuss the rel- 
evant physics of the two-impurity model. The study of 
the critical point for two impurities interacting with the 
Hamiltonian Eq.(^) starts with a reduction of a three- 
dimensional problem to a one-dimensional one P,p|,|l0| . 
This is done by introducing orthogonal II? fermionic op- 
erators, 



f/> lt a(k) = I dn— 



1 



N e (k) 



cos(k ■ R/2) ± 



N (k) 



sin(k ■ R/2) 



o(k), 



(11) 



where we have integrated over spherical angle, and Then the two-impurity model can be rewritten in ID 

I : — form. While the transformation itself is exact, the k- 

, \ / sin(kR) , x 

N e , (k) = \/l±—^. (12) 
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dependent couplings arise. Then these couplings are as- 
sumed to be /c-independent for excitations in the vicinity 
of the Fermi surface, an assumption questioned in Ref. 

. The main drawback of such assumption is that the 
RKKY exchange interaction then does not have a cor- 
rect oscillating character. The reason is that the main 
contribution to RKKY interaction comes from energies 

H int = iA$(0)<r-Si&(0) + i 



~ vp/R, where vp is the Fermi velocity, R is the dis- 
tance between impurities. Apart from this, however, this 
procedure is fully justified. One can view the process 
of cutoff renormalization in two stages. During the first 
stage, the cutoff A is reduced below vp/R, A -c Vp/R, 
and oscillating RKKY interaction appears as a result: 

2(0)<r-S 2 &(0)+ISi-Sa. (13) 



Here A = Jp is the Kondo coupling constant and 
I = V cos(2kpR)/R 3 is the usual RKKY term, V = 
pj 2 /(167r). Note that both impurities are now at the 
same place in the ID model. During the second stage 
the ID Hamiltonian Eq.(|l|) is renormalized. The wave 
functions tpi and ip2 are not orthogonal; they can be writ- 
ten in terms of orthogonal ipi and tp2 as follows: 



where 







ip 2 = Wi + aip 2 , 



N e {k F ) + N (k F ) 

2 ' 
N e (k F ) - N {k F ) ^ sin(kpR) 

2 ~ 2k F R 



(14) 
(15) 



(16) 
(17) 



This orthogonal basis differs from the usual "even" and 
"odd" parity wavefunctions. It is seen explicitly in this 
basis that for kpR 3> 1 channel mixing terms play mi- 
nor role and can be neglected. Other marginal operators 
appear in the interaction Hamiltonian Eq.(13). Those, 
however, are suppressed by factors of kpR, and therefore 
can be neglected. 



V. HIGH-TEMPERATURE RESULTS. 



We can now analyse the nf corrections in a Kondo alloy 



of spin-1/2 impurities using the II? Hamiltonian Eq.(13). 
At high temperatures (T ^> T K ) this can be done by in- 
cluding the RKKY term in Eq. ([l3]) in the non- interacting 
Hamiltonian, and perturbing in the Kondo interaction. 



Instead of averaging over distances as in Eq.(lO) one can 
do averaging over the strength of the RKKY coupling 
p2[ . Indeed, the interaction I(R) is a product of the 
smooth function R~ 3 times the rapidly oscillating func- 
tion cos(2kpR). It is convenient to carry out the integra- 
tion in formula Eq.([lo|) by first averaging the integrand 
over the period of oscillations. Changing the order of 
averaging and integration, we obtain: 



i?(2) 



3T 



v.p.J ^(F 2i (y,X)-F, 



2,(0, A)), (18) 



where y = I(R)/T, v. p. stands for "principle value" . The 
behavior of the bulk spin susceptibility and the heat ca- 
pacity of a Kondo alloy at temperatures T 3> Tk is well 
known: 



X(T) = — (1 - A e// (T)) 

C(T) = -y ni \ e}f {T)\ 



(19) 



where A e //(T) ~ A + A 2 In [^] is the effective coupling 
constant at energy scale T. To the lowest orders in the 
Kondo coupling we get the following result for the nf 
correction to the heat capacity: 



SC 



3 T 



1 - A In 



while for the magnetic susceptibility we find: 



S\ = k (5Mb) 



2 n i ^0 

3T 2 



1 - A - 2A In 



where 



v.p. 



dyl 



-oo y 



-0.501447, 



(20) 



(21) 



(22) 



and A < vp/R ~ v F (T K /V ) 1/3 is the cutoff of the ID 
model. Note that the small-bandwidth effective theory of 
the two-impurity Kondo model is, of course, only valid for 
calculating low energy properties. Thus the above results 
are valid in the range T K < T <S A. Since the RKKY 
interaction is relevant at weak coupling, it gives rise to 
power law singularities in thermodynamic functions at 
high temperatures fl2|| . The marginal Kondo interaction 
gives rise to weaker logarithmic terms in Eqs. (pp|) , (^T|) . 
As we will see below in the next section, these singulari- 
ties are cut off at temperatures of the order of Tk- The 
parameter of the virial expansion is UiVo/T <C 1. The 
major contribution to the nf term comes from distances 
R - (y /T) 1/3 , where the amplitude of the RKKY term 
is of the same order as temperature. The virial expansion 
fails at high enough densities when the average distance 
between impurities becomes comparable to R. Note that 
the second virial coefficients for the spin susceptibility 
and the heat capacity Eqs.(pp|),([2l|) cannot be written 
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simply in terms of an effective coupling. Averaging over 
distances in the virial expansion correspond to different 
values of RKKY exchange I, which can be quite large, 
ft is clear that I ruins the simple dependence of phys- 
ical quantities on A and bare coupling. Consider what 
happens when we take I either very large or very small, 
continuing to treat I exactly and doing perturbation the- 
ory in the Kondo couplings. If / is very small (compared 
to T) we should be able to ignore it and then we just 
get two decoupled single impurity Kondo problems. In 
this limit we certainly recover the standard reduction of 
A-dependence to renormalized coupling. On the other 
hand, if T <C |/| (/ < 0) then, at low energies, we have 
essentially a two-channel S = 1 Kondo problem. Again 
we expect a standard reduction of the dependence on 
A corresponding to the beta function for that problem. 
However, the evolution of effective Kondo coupling with 
energy scale is different in these two limits (two decou- 
pled one-channel S — 1/2 Kondo problems versus one 
two-channel S=l Kondo problem). Thus, as we lower our 
cut-off past / the nature of the dependence on the cut- 
off changes. This is to be expected on general grounds. 
When we have a physical energy scale, like / in the prob- 
lem we don't in general get simple dependence on the 
cut-off scale. In a massive quantum field theory effective 
coupling constants are only useful at energy scales large 
compared to the mass. 



VI. CORRECTIONS AT LOW TEMPERATURES. 

Let us now consider this problem at low temperatures. 
Averaging over the distances between two impurities in- 
volves integration over all RKKY couplings (see Eq.([l8|)), 
which makes analytical analysis extremely difficult. In- 
stead let us consider first the SU(N) generalization of 
the two-impurity model, for which explicit results for the 
linear term in the specific heat (7) and the bulk magnetic 
susceptibility were obtained as a function of RKKY 
coupling strength. For a U — 00 SU (N) generalized two- 
impurity Anderson model with RKKY interaction they 
find, in the large- A limit: 



7 



Ntt 2e- BS cosS 
3 {l-A 2 )T i 



(23) 



K 



where A ~ sinikpR) / kpR, B ~ cos(kFR)/kpR, S is re- 
lated to the scattering phase shift in even and odd chan- 
nels, S e . — tt/2 ± S, N is the SU(N) expansion param- 
eter, T K is the Kondo temperature. S is found from the 
solution of appropriate mean field equations. To find the 
nf correction to 7 and x, we have to average Eq.(p3|) over 
distances, or over RKKY interaction. Since for the aver- 
age distance UfR 1, B — > 0, A — > 0. For these values 
of parameters the model has a first-order phase transi- 
tion at (I/T K ) C = 8/tt, with 5-0 when (I/T K ) < 8/ir 



and 5 ~ 7r/2 if (I/Tk) > 8/tt. We can apply the same 
averaging technique as in Eq.dia) (with y — I(R)/Tk)- 



r 



^v. P .[ d 4( X2i (y,X)-X 2i (0,X)). (24) 
3T K J y 2 



Performing the integration, we easily find: 



7 = 7T X = Ni 



7TA 

37V 



1 



(25) 



where the first term in Eq.(p5|) corresponds to the or- 
dinary Kondo effect in the SU{N) model. One impor- 
tant result of the virial expansion is to find the density 
at which the single-impurity approximation no longer 
works. As one can see from Eq.(|25]), this happens at 
fii ~ Tk/Vq, or when the Kondo temperature becomes 
comparable to the amplitude of the RKKY interaction 
at the average inter-impurity distance. 

The mean field N = 00 solution suffers from a number 
of defficiencies. In particular, it was shown [|6| JTl1Jl0tl that 
for two spin- 1/2 impurities the quantum phase transition 
is second, not first order. The critical point for two spin- 
1 /2 impurities is only present if the model is particle- hole- 
symmetric. Under these conditions it was found [pd|Jl0[ 
that the specific heat coefficient 7 diverges at the critical 
point as \I/T K - (I/T K ) C \- 2 , {I/T K ) C - 2.2, while the 
uniform spin susceptibility X does not. Let us find out 
whether this divergence produces any singularity in the 
thermodynamic quantities. Unlike in the N — 00 model, 
where there is always a first order phase transition, One 
can identify possible singularities at y ~ and y ~ y c 
from Eq.(|24|). At y = I/Tk —> the two Kondo impuri- 
ties are independently screened. Since impurity spin at 
the Kondo fixed point should be replaced by the local 
spin density, 



Si cx 



ufV>1(0Wl(0) 
T K 



(26) 



so that the RKKY interaction can be written as follows: 

Hrkky « -^(V>i(0Wi(0) • </4(0W 2 (0)). (27) 

Perturbing in Hrkky, on e easily finds that the singu- 
larity in 7 and x at ?/ — > is logarithmic, i.e. principle 
value-integrable. For y ~ y c the bulk susceptibility x 
does not have a singularity, whereas 7 has a power law 
divergence ~ {y — y c ) 2 - The form of this singularity in 
the vicinity of y c was investigated ealier using bosonisa- 
tion by Gan p0| , who showed that it is effectively cut 
off at finite temperatures at Sy cx \/T. This produces a 
square root singularity in 7, 



Nt n z V 



T K ^/TTk 



(28) 



Since x does not diverge at the critical point, such singu- 
larity is absent for the bulk susceptibility. It is known jjj] 
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that in the absence of particle-hole symmetry there are 
relevant perturbations near the zero-temperature criti- 
cal point of the two-impurity Kondo model. Such per- 
turbations tend to wipe out the singularity at low tem- 
peratures, so it is not clear whether this effect could be 
observed in real materials. 



VII. CONCLUSIONS. 
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To summarize, we have considered cluster expansion 
in "dense" Kondo alloys. This expansion ties the two- 
impurity model, which has been investigated in detail, 
to experiment. At high temperatures nf corrections pro- 
duce additional 1 /T and logarithmic singularities in ther- 
modynamic quantities. At low temperatures the n 2 term 
becomes large when the RKKY interaction at the av- 
erage distance between impurities becomes of the same 
order as the Kondo temperature. Since the parameter 
of expansion in density is iii/n c , where n c — Tk/Vo oc 
A -2 exp(— 1/A) (A = pj is the dimensionless Kondo cou- 
pling constant) is exponentially small, deviation from the 
single-impurity behavior could be observed in dilute al- 
loys with magnetic impurities. We have also seen that 
the presence of the intermediate fixed point in the two- 
impurity problem could lead to a mild low-temperature 
singularity in the specific heat, but not the bulk mag- 
netic susceptibility. Further numerical investigation of 
this problem is desirable. In some rare earth alloys, such 
as Yb x Yi_ x CuAl |EJ, the behavior of thermodynamic 
functions is linear in x up to very large x ~ 0.9. This ob- 
servation could be related to the fact that these ions are 
often described by the Coqblin- Schrieffer model (N=8 in 
case of Yb 3+ ), where the RKKY interaction is suppressed 
(0(1/Y 2 ) @). 

Our calculation of the second virial coefficient shows 
that the density of impurities at which multi-impurity 
effects become important is higher than that given by 
Nozieres' exhaustion principle. According to this prin- 
ciple, at these densities (n < n c ) there are not enough 
states to screen all impurities. Since Nozieres' principle 
is multi-impurity in nature, it may be necessary to study 
the convergence properties of the expansion itself rather 
than the behavior of a given term to claim its failure. 

Finally, we note that the reason why the large Kondo 
length scale is absent in the problem is purely geomet- 
ric. The channel mixing terms between two impurities 
are large in the one-dimensional problem. However, in 
3D they become suppressed by a small factor of I/UfR, 
and the physics of the two-impurity model (and related 
nf term in the cluster expansion) is fully determined by 
the competition of the RKKY and Kondo effects. 
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